Geometry in transition: A model of emergent geometry 
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We study a three matrix model with global SO(3) symmetry containing at most quartic powers 
of the matrices. We find an exotic line of discontinuous transitions with a jump in the entropy, 
characteristic of a 1st order transition, yet with divergent critical fluctuations and a divergent 
specific heat with critical exponent a = 1/2. The low temperature phase is a geometrical one with 
gauge fields fluctuating on a round sphere. As the temperature increased the sphere evaporates in 
a transition to a pure matrix phase with no background geometrical structure. Both the geometry 
and gauge fields are determined dynamically. It is not difficult to invent higher dimensional models 
with essentially similar phenomenology. The model presents an appealing picture of a geometrical 
phase emerging as the system cools and suggests a scenario for the emergence of geometry in the 
early universe. 
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i— h ' Our understanding of the fundamental laws of physics has evolved to a very geometrical one. However, we still 
I ^ . have very little insight into the origins of geometry itself. This situation has been undergoing a significant evolution 
in recent years and it now seems possible to understand classical geometry as an emergent concept. The notion of 
geometry as an emergent concept is not new, see for example [l[ for an inspiring discussion and 0, Q for some recent 
i— i' ideas. We examine such a phenomenon in the context of noncommutative geometry Q emerging from matrix models, 
by studying a surprisingly rich three matrix model 0, 0, 0] . The matrix geometry that emerges here has received 
attention as an alternative setting for the regularization of field theories [1, 0, HI HI and as the configurations of DO 
branes in string theory [l2|, EH- Here, however, the geometry emerges as the system cools, much as a Bose condensate 
or superfluid emerges as a collective phenomenon at low temperatures. There is no background geometry in the high 
temperature phase. The simplicity of the model, in this study, allows for a detailed examination of such an exotic 
transition. We suspect the asymmetrical nature of the transition may be generic to this phenomenon. 

We consider the most general single trace Euclidean action (or energy) functional for a three matrix model invariant 
under global SO(3) transformations containing no higher than the fourth power of the matrices. This model is 
t-H , surprisingly rich and in the infinite matrix limit can exhibit many phases as the parameters are tuned. We find that 
generically the model has two clearly distinct phases, one geometrical the other a matrix phase. Small fluctuations 
in the geometrical phase are those of a Yang-Mills theory and a scalar field around a ground state corresponding to 
a round two-sphere. In the matrix phase there is no background spacetime geometry and the fluctuations are those 
of the matrix entries around zero. In this note we focus on the subset of parameter space where, in the large matrix 
limit, the gauge group is Abelian. 

For finite but large N, at low temperature, the model exhibits fluctuations around a fuzzy sphere In the infinite 
N limit the macroscopic geometry becomes classical. As the temperature is increased it undergoes a transition where 
the entropy jumps, yet the model has critical fluctuations and a divergent specific heat. As this critical coupling is 
approached the fuzzy sphere radius expands to a critical radius and the sphere evaporates. The neighbourhood of the 
critical point exhibits all the standard symptoms of a continuous 2nd order transition, such as large scale fluctuations, 
critical slowing down (of the Monte Carlo routine) and is characterized by a specific heat exponent which we argue 
is a = 1/2, a value consistent with our numerical simulations. In the high temperature (strong coupling) phase the 
model is essentially a zero dimensional Yang-Mills theory. 

By studying the eigenvalues of operators in the theory we establish that, in the matrix phase, the matrices D a are 
characterised by continuous eigenvalue distributions which undergo a transition to a point spectrum characteristic of 
the fuzzy sphere phase. The point spectrum is consistent with D a — L a /R where L a are su{2) angular momentum 
generators in the irreducible representation given by the matrix size and R is the radius of the fuzzy sphere. The full 
model received an initial study in Q while a simpler version invariant under D a — ► D a + Al arises naturally as the 
configuration of DO branes in the large k limit of a boundary Wess-Zumino-Novikov-Witten model and has been 
studied numerically in The interpretation here is novel, as are the results on the entropy and critical behaviour 
and the extension to the full model. 
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Let D a , a 



1,2,3, be three NxN Hermitian matrices and let us consider the action 
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so that a 4 = (3 plays the role of inverse temperature and we fix c 2 = J2 a ^a — (^ 2 — m this study [23J]. The 
absolute minimum of the action is given by D a = L a . Expanding around this configuration, with D a — L a + A a , 
yields a noncommutative Yang-Mills action with field strength 



F ab = i[L a , A b ] - i[Lb, A a ] + e a b c A c + i[A a , A b 
and gauge coupling g 2 — 1/a 4 . As written the gauge field includes a scalar field, 

1 A 2 

$ = - {x a A a + A a X a + —=) , 

2 Jci 



(2) 



(3) 



as the component of the gauge field normal to the sphere when viewed as embedded in R 3 where x a — L a j \fc2. 

In the large N limit, taken with a and m held fixed the action for small fluctuations becomes that of a U(l) gauge 
field coupled to a scalar field defined on a background commutative two-sphere Iff. For large m the scalar field is not 
excited. The model with m — appears as a low energy limit of string theory [131 ] . 

One can see the background geometry as that of a fuzzy sphere [l4| by noting that the x a satisfy 
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and the algebra generated by products of the x a is the algebra of all N x N matrices with complex coefficients. The 
geometry enters through the Laplacian [§] 



£ 2 - = [L a ,[L a ,-]] , 
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which has the same spectrum as the round Laplacian on the commutative sphere, but cut off at a maximum angular 
momentum L = N — 1. The fluctuations of the scalar have this Laplacian as kinetic term. 

The ground state is found by considering the configuration D a = 4>L a where ip plays the role of the inverse radius 
of the sphere. The effective potential [^.Il5| obtained by integrating out fluctuations around this background is given, 
in the large TV limit, by 
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It is unbounded from below at the origin, but, for sufficiently large a, exhibits a local minimum near </> = 1 which 
disappears for a < a* . The critical curve <S* is determined from the point at which the real roots of dV c s /dip = 
merge and disappear. This occurs at 
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and gives the critical curve 
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where t = m 2 (l 



?). This expression interpolates between a* — (8/3) 3 ^ 4 with 0„ = 3/4 for m = and a* = 

(8/to 2 ) 17 *, with (f)* = 1/V2 for large m. Thus, as the system is heated, the radius, R — <j)~ 1 , expands form R = 1, at 
large a to R^ 1 = (f)* at a*. When a < a* the fuzzy sphere solution no longer exists, and, in effect, the fuzzy sphere 
evaporates as the radius approaches i?*. 
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Furthermore, defining S =< S > /N 2 , we have 
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which for the fuzzy sphere phase becomes 
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with <fi given by the local minimum of V e ff . Expanding in the neighbourhood of the critical point gives 




1(02+2^) 



(/3-/3 c ) 1/2 + ... 



(11) 



and predicts that the transition has a divergent specific heat with exponent a — 1/2. For m = we find S* = 5/12 
while in the matrix phase it takes the constant value S = 3/4 hence AS — 1/3, i.e. the entropy has a discrete 
jump of 1/9 per degree of freedom (there are 3 matrices) across the transition. These fractions are in excellent 
agreement with the values obtained numerically suggesting that a more complete exact solution to the model, at least 
with m = 0, may exist. 

Monte Carlo Simulations: The model ([I} with m = was studied in [5[ and an initial study of the full model 
was reported in 0] • In Monte Carlo simulations we use the Metropolis algorithm and the action ([I} . The errors were 
estimated using the jackknife method. 

We observe that S collapses well for different N. Defining a s as the value of a at which curves of the average value 
of the action, < S >, for different N cross, gives a good estimate of the location of the transition. 

For m — we observe (figure [J) a divergence in the specific heat, C v :—< (S— < S >) 2 > /N 2 , as the critical 
value a s is approached from above, i.e. from the low temperature phase. For a > a s the model is in the fuzzy sphere 
phase and C v rapidly approaches C v = 1 as a is increased. For a < a s the model is in a matrix phase and C v = 0.75, 
retaining this constant value right from the transition at a s to a = 0. For large N we observe that the location a max 
of the peak in C v and the minimum a m in coincide to the numerical accuracy we have explored and agree well with 
the coupling, a s , at which < S > intersect for different values of N. The critical coupling determined either as a max 
or a m i n from C v , or as a s gives good agreement with ((5]) and the numerical value found in || . 

The behaviour when m ^ is qualitatively similar to m — 0. We observe again that, for large a, the model is in 
the fuzzy sphere phase where C v = 1. 

As a is decreased for fixed TV, C v goes through a peak at a max , with the peak value decreasing as m 2 is increased. 
As a is further decreased from a max , C v goes through a minimum at a m i n and then increases again, while asymptoting 
to the value 0.75 at a = 0. As the value of m 2 is increased, our numerical study confirms that the fuzzy sphere-matrix 
model transition is shifted to lower values of a. Extrapolating m 2 — > oo we infer that the critical coupling goes to 
zero. This limiting model should be related to the fuzzy Yang-Mills model without scalar field [3. By extrapolating 
the measured values of a max and a m i n to N = oo wc obtain the critical value a c . This agrees with a s and hence we 
infer both are detecting the matrix-to-sphere phase transition. However, a c deviates measurably from ([5]); e.g., for 
rn 2 = 200 a c = 0.403 ± 0.014 while © gives &* = 0.446982. Thus, though © gives a good indicative value for the 
transition it is not precise. Our results are summarised in a phase diagram in figure O 

We further observe that the jump in entropy decreases as m is increased and the nature of the transition seems 
to change. Our theoretical analysis (fTTj) still indicates a divergent specific heat with exponent a = 1/2 but with a 
narrowing critical regime as m is increased. However, we find no direct evidence for such a regime in the numerical 
data for large m. We have not been able to determine with any precision where the transition becomes continuous; 
however, the entropy jump seems to disappear at m 2 ~ 20. Also the data for a s and a c separate in this region (see 
figure[2|) indicating a possible multi-critical point. It is also possible that beyond this point the transition is in fact 3rd 
order, a behaviour typical of many matrix models; however, the persistence of the critical line as determined by a s 
the crossing point of the < S > curves suggests the transition is 2nd order continuous, consistent with the theoretical 
analysis. We know of no other model that exhibits transitions of the type presented here and further numerical and 
theoretical study is warranted. 

More detail on the structure of the phases can be obtained from the distribution of eigenvalues of observables. 
Here we focus on D$ (by rotational symmetry all matrices have the same eigenvalue distribution) . The characteristic 
behaviour of the distributions of eigenvalues in the fuzzy sphere and matrix phases is illustrated in figures [3] and 0] 
respectively. Thus we see that, as one crosses the critical curve in figure [21 the eigenvalue distribution of D% undergoes 
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FIG. 1: The specific heat for the model So as a function of coupling or inverse temperature (3 = a 4 . The theoretical curve is 
obtained from the local minimum of the effective potential Eq. @ 
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FIG. 2: The phase diagram. The theoretical curve is given by Eq. Q. The specific heat exhibits a divergence for small m, 
with the transition becoming smooth around lnm 2 ~ 2 where the curves a s and a c separate. 
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a transition from a point spectrum to a continuous distribution symmetric around zero, the latter can be fit by the 
one-cut distribution [20| for a single matrix in a quartic potential. 
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FIG. 3: The eigenvalue distribution for D3, with TV = 24, in the fuzzy sphere for a — 5 and m 2 
D 3 = L 3 . 
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Conclusions: We have found an exotic transition in a simple three matrix model. The nature of the transition 
is very different if approached from the high or low temperature side. The high temperature phase is described 
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by Yang-Mills in zero dimensions so there is no background spacetime geometry. The system in this phase is well 
approximated by 3 decoupled matrices each with the same quartic potential. The value C' v — 0.75 coincides with the 
specific heat of 3 independent matrix models with quartic potential in the high temperature limit and is therefore 
consistent with this interpretation. 

As the system cools a geometrical phase condenses and at sufficiently low temperatures the system is described by 
small fluctuations on a two sphere of a U(l) gauge field coupled to a massive scalar field. The critical temperature is 
pushed upwards as the scalar field mass is increased (see figure [2]). We also find evidence that the UV/IR mixing [22| 
typical of many non-commutative field theories is linked to the geometrical instability and disappears as the scalar 
mass is increased. Once the geometrical phase is well established the specific heat takes the value 1 with the gauge 
and scalar fields each contributing 1/2. 

The model of emergent geometry described here, though reminiscent of the random matrix approach to two dimen- 
sional gravity Q is in fact very different. The manner in which spacetime emerges is also different from that envisaged 
in string pictures where continuous eigenvalue distributions j2flor a Liouville mode 21 1 give rise to extra dimensions. 
It is closely connected to the DO brane scenario described i n 1 1 21] and our model with m = is a dimensionally reduced 
version of a boundary WZNW models in the large k limit [13j . It is not difficult to invent higher dimensional models 
with essentially similar phenomenology to that presented here (see 17 1, [l9| and [Hj]). 

We believe that the model presented above gives an appealing picture of how a geometrical phase might emerge as 
the system cools and suggests a very novel scenario for the emergence of geometry in the early universe. In such a 
scenario the temperature can be viewed as an effect of other degrees of freedom present in a more realistic model but 
not directly participating in the transition we describe. 
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